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Abstract —Visible light communication (VLC) using light- 
emitting-diodes (LEDs) has been a popular research area re¬ 
cently. VLC can provide a practical solution for indoor position¬ 
ing. In this paper, the impact of multipath reflections on indoor 
VLC positioning is investigated, considering a complex indoor 
environment with walls, floor and ceiling. For the proposed 
positioning system, an LED bulb is the transmitter and a photo¬ 
diode (PD) is the receiver to detect received signal strength (RSS) 
information. Combined deterministic and modified Monte Carlo 
(CDMMC) method is applied to compute the impulse response 
of the optical channel. Since power attenuation is applied to 
calculate the distance between the transmitter and receiver, the 
received power from each reflection order is analyzed. Finally, 
the positioning errors are estimated for all the locations over 
the room and compared with the previous works where no 
reflections considered. Three calibration approaches are proposed 
to decrease the effect of multipath reflections. 

Index Terms —Indoor positioning, visible light communication, 
multipath reflections, impulse response, received signal strength. 


I. Introduction 

LOCATION BASED SERVICES (LBS) have become a 
popular research topic for several years which provide 
users with current locations and related services. For outdoor 
environment. Global Positioning System (GPS) provides satis¬ 
factory services such as localization, navigation and displaying 
surrounding traffic conditions. For indoor environments, GPS 
technology is not applicable since a satellite signal suffers 
from severe attenuation when passing though solid walls. In 
recent decades, several methods have been proposed to realize 
indoor positioning with the help of technologies such as ultra¬ 
wide band (UWB), wireless local area network (WLAN), 
Radio-frequency identification (RFID), Bluetooth and cellular 
system m. 

Light-emitting-diode (LED) technology has been develop¬ 
ing very rapidly in recent decades. It not only provides people 
with economical and efficient illumination and a long service 
time, but also paves the way for smart lighting and visible 
light communication (VLC) 0-01. As a strong candidate for 
high-speed wireless networks of the next generation, visible 
light communication exhibits many advantages over conven¬ 
tional radio frequency (RF) communication. First, visible light, 
together with infrared and ultraviolet spectral band, provides 
unregulated and unlimited bandwidth as a practical solution to 
the current spectrum crunch issue 0, 0. Second, considering 
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that light waves are unable to penetrate through solid walls 
and are confined in an individual room, the band reuse among 
different rooms is accomplished and physical layer security 
for the communication system is guaranteed. Third, VLC can 
be widely applied in many RF sensitive environments such 
as mines, power plants and hospitals due to the fact that light 
waves never generate any electromagnetic interference. Fourth, 
so long as illumination infrastructure exists, VLC is applicable 
so that the hardware cost is decreased. 

The application of VLC technology for the indoor position¬ 
ing has been extensively studied as an available solution for 
the LBS I9l-fl7l. LED light sources act as transmitter and 
receiver is a photo-diode (PD) or an image sensor collocated 
with a user. Several approaches have been proposed to realize 
visible light positioning. In one approach, an image sensor is 
used to obtain angulation algorithm to calculate the receiver 
position based on angle-of-arrival (AOA) information and 
rotation matrix |[T0|. In this method, colored LEDs are used 
to help the image sensor to distinguish between different light 
sources. Scene analysis is another approach to obtain the 
receiver position. Features of each location are collected as the 
fingerprints in the offline stage. In the online stage, the features 
of current location are measured and by matching those with 
offline fingerprints, location of receiver is estimated ITU . In 
this paper, we employ a commonly used algorithm where the 
RSS information is first detected by a PD, and then distance 
between transmitter and receiver is calculated. The lateration 
algorithm is finally applied to estimate receiver coordinates 

m. 

In addition to the above methods, other technologies are 
introduced in VLC system to improve the positioning per¬ 
formance. Zigbee technology can be combined with VLC 
to realize long distance positioning Gl- In OH, with the 
assistance of a 6-axes sensor (geomagnetic sensor and gravity 
acceleration sensor), a switching estimated receiver positioning 
system is proposed to achieve higher accuracy. Hybridizing ac¬ 
celerometer was proposed in ED to realize three dimensional 
positioning without knowledge of receiver height. Gaussian 
mixture sigma point particle filter can be further employed 
to achieve high positioning accuracy and prevention of large 
deviations 110. 

In the literature, line-of-sight (LOS) channels have been 
considered without taking account of multipath reflections 
in analysis of positioning performance. However, transmitted 
signal introduces multipath reflections as it bounces off walls, 
ceiling and floor where the transmitter is a wide-beam LED 
source, and the receiver having a finite field-of-view (FOV) 
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captures reflected photons from room surfaces. In this pa¬ 
per, the effect of multipath-induced distortion on positioning 
accuracy of indoor VLC positioning systems is investigated. 
Several methods have been proposed to approximate impulse 
response of an indoor optical wireless channel. In lUSl . Barry 
et al. proposed a deterministic algorithm that partitions a room 
into many elementary reflectors and sums up the impulse 
response contributions from different orders of bounces. As 
this method is recursive, it takes a significant amount of 
computation time. Monte Carlo ray tracing approach is an 
alternative way to calculate the impulse response where rays 
with identical optical power are traced from the source ED- 
The direction of rays are generated by a probability density 
function (PDF) modeled using a Lambertian pattern. When 
the rays hit reflecting surfaces, new rays with reduced power 
are generated from the impact point with the same PDF. This 
method suffers from the fact that it needs a very large number 
of rays while only a small portion of rays will finally reach 
the PD. To alleviate this issue, modified Monte Carlo (MMC) 
approach is proposed which exploits each ray several times 
instead of only once f20ll . Although this method is fast, it 
introduces some variance due to random direction of the rays. 
In El. Alqudah and Kavehrad proposed a new approach 
to characterize diffuse links in a multiple-input multiple- 
output (MIMO) system. In this paper, we use the methodology 
developed recently in lf22l . referred as combined deterministic 
and modified Monte Carlo (CDMMC), taking advantage of 
both methods to simulate impulse responses of indoor optical 
wireless channels. 

The rest of the paper is organized into five sections. In 
Section II, the system model and CDMMC approach used to 
calculate the impulse response are briefly discussed. In Section 
III, the positioning algorithm is described to estimate receiver 
coordinates. In Section IV, the effect of multipath reflections 
on the positioning system is thoroughly investigated. In Sec¬ 
tion V, three calibration methods are applied to increase the 
positioning accuracy. Section VI finally concludes the paper. 

II. Multipath Analysis 
A. System Design 

A typical indoor optical positioning system is shown in Fig. 
1 where the LED bulbs as transmitters are installed on the 
ceiling of the room, and the receiver is located 1.2 m from the 
floor. Each of the transmitters has an identification (ID) code 
denoting its coordinates. On-Off-keying (OOK) modulation is 
used to modulate the LED bulbs. Six reflection surfaces of the 
room, i.e., four walls, one ceiling and one floor are assumed 
to be perpendicular with one another. 

As shown in Table 1, the reflection coefficients are assumed 
to be fixed considering the material of the room surface and 
the 420 nm light source. The transmitters in the algorithm are 
treated as point sources and located at the height of 3.3 m 
considering practical installations. The distance between the 
transmitters is 2 m. As the transmitters are facing downwards, 
the azimuth angle is 0° and the elevation angle is -90°. The 
receiver is facing upwards, thus the azimuth angle is 0° and 
the elevation angle is +90°. The receiving area of the PD is 



TABLE I 

Room Configuration 


Room dimensions 

Reflection coefficients 

length: 6 m 
width: 6 m 
height: 3.5 m 

Pwall • 0-66 
PCeiling • 0-35 
PFloor ■ 0-60 

Horizontal coordinates of LED bulbs: 

(1.0 , 1.0) (1.0 , 3.0) (1.0 , 5.0) (1.0 , 7.0) 

(3.0 , 1.0) (3.0 , 3.0) (3.0 , 5.0) (3.0 , 7.0) 

(5.0 , 1.0) (5.0 , 3.0) (5.0 , 5.0) (5.0 , 7.0) 

(7.0 , 1.0) (7.0 , 3.0) (7.0 , 5.0) (7.0,7.0) 


1x10 4 m 2 , with 70° FOV. These parameters are summarized 
in Table 1. 


B. Channel Access Method 

As the LED bulbs transmit their coordinates information 
independently, the signals cannot be retrieved if they interfere 
with each other at the receiver. To solve this problem, time 
division multiplexing (TDM) is applied as the channel access 
method (23). All the transmitters have synchronized frames 
and occupy different time slots in one frame period to send 
their signals. The frame structure is presented in Fig. 2. When 
one LED bulb transmits the ID information, other LED bulbs 
emit constant light intensity for the illumination purpose. 


LED 1 

ID code 

constant intensity 






LED 2 

constant 

intensity 

ID code 
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• 
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LED N 

constant intensity 

ID code 
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Time Time Time 

Slot 1 Slot 2 Slot N 


Fig. 2. Frame structure of the positioning system for one period. 
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C. Impulse Response Analysis Method 


The CDMMC method combines the deterministic and the 
MMC methods to take advantage of both. The deterministic 
method is employed to calculate the contribution of the first 
reflections for high accuracy ED, while the second and higher- 
order reflections are calculated by the MMC method to achieve 
a high computational speed m- 

The first step is to divide the room surfaces into many small 
square elements each of which has an area of 1 x 10” 4 m 2 , 
equal to the PDs receiving area. 

Second, the PD together with these elements act as the 
receivers. The received power is obtained as 

Preceived (0) = J ff(0)P SO urce ( ° ) , d) 


where P sOU rce < ' 0 ' ) is the power emitted from each LED bulb. In 
([]}, H (0) is the channel DC gain El . and P SO urce^ is the 
received power for the LOS link. 

Third, each of these small elements is considered as a point 
source whose power is calculated as 


p (t) 
1 source 


— -^received 


Psurface i 


( 2 ) 


where ^surface is the reflection coefficient of ceiling, floor or 
walls listed in Table I. 

In the fourth step, the received power of the first reflections 
is calculated by considering the small elements and PD as the 
receivers again which can be expressed as 

Prece,ved (1) =P(0)P SOU rce (1) . (3) 


The MMC method is then employed such that 10 random rays 
are generated from each small element where the PDF of the 
rays’ directions follows 

777 -L 1 

f(a : /3) = ^-cos m (a). (4) 

In (0]), a is the angle between z-axis and the ray vector shown 
in Fig. 3, /3 is the angle between projection of the ray vector 
on the X — Y plane and x-axis, and m is the Lambertion order. 
In Fig. 3, the origin point is each ray’s point source and the 
X — Y plane represents the surface plane of the source. Note 
that (0]) is independent of 3 .These rays arrive at the surface of 
the room and then the received power of secondary reflections 
is calculated as 


Psource ' — -^received ' Psurface/10. 


(5) 






LOS 


First 

Reflections 


Second 

Reflections 


Subsequent 

reflections 




• • • 


Fig. 4. Flow diagram of CDMMC method. 


(i) 

logic 0. The relation between P r and I\ can be expressed 
as El 


The impact points of these rays on the room surfaces are 
considered as the new sources, and then new rays are generated 
to calculate the subsequent reflections, recursively. Fig. 4 
shows the flow diagram of the CDMMC method. The impulse 
response of the channel is computed by adding up all the 
contributions from LOS and those of reflections. 

III. POSITIONING ALGORITHM 
A. Distance Estimation 

After the impulse response is computed, received signal 
strength for transmitter i is obtained as P r ^, we assume that 
the transmitted power, i.e. P t , is 5 W for logic 1 and 3 W for 


P r W = ™ + .l Acos m (cp) T s (ip) g (ip) cos (ip) P t , (6) 

ZTTdi 

where <p is the irradiance angle from the transmitter to the re¬ 
ceiver, ip is the incident angel, and T s (ip) is the transmittance 
function of the optical filter. In (H, g(ip) is the compound 
parabolic concentrator gain of the receiver El . and A is the 
detector physical area. Table II presents the numerical value of 
these parameters assumed throughout this paper. Considering 
transmitters and receiver’s orientation, cos (</>) = cos (ip) = 
h/di where h is the vertical distance between transmitter and 
receiver. Therefore, the distance di between the transmitter i 
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TABLE II 

Parameters of the transmitter and receiver 


Transmitters (Sources) 

Receiver 

Wavelength: 420 nm 
Height ( H ): 3.3 m 
Lambertian mode (m): 1 
Elevation: -90° 
Azimuth: 0° 

Area (A): lxl0 - ' 1 m /: 
Height ( h ): 1.2 m 
Elevation: +90° 
Azimuth: 0° 

FOV (® c ): 70° 
T.W):l 

SW): 2.5481 


and the receiver is calculated as 


di = 


(m + 1) AT S {ip) g {ip) P t h 2 


2t tP 


(i) 


(7) 



The horizontal distance is finally estimated as 


Fig. 5. Birds-eye view of the system model. 


= \ d 2 — h 2 = 


\ 


' (m + 1) AT S {ip) g {ip) P t h 2 


2t rP 


M 


-h 2 . 

( 8 ) 


is obtained and by setting ^ to zero, the following equation 
is derived 


— 2A t B + 2A t AX = 0. (15) 


B. Least Square Estimation 

The receiver coordinates {x, y) correspond to the transmitter 
coordinates {xi, yP) based on the following group of equations: 

' {x - zi) 2 + {y- yif = r\ 

{x - x 2 f + {y- y 2 ) 2 = r\ 

< . , (9) 

\ x ~ Xnf + {y- ynf = r 2 

where n is the number of transmitters whose signal can be 
detected by the receiver. By subtracting the last three equations 
from the first one, we obtain 


r {xi -x 2 )x + {yi -y 2 )y 
= ( r l -r{-xl+x\-yl + y\) /2 
< : 

(zi -x n )x + {y 1 -y n )y 
.= ( r l - r l - X l + X \ ~ Vi + Vl) / 2 

which can be expressed in the matrix format as 

AX = B 

where 


A = 


x 2 - xi y 2 y i 


B = - 

2 


_x n - xi y n - y i 
(r 2 - rl) + (xl + yl) - (x\ + yf) 

(r\ - rl) + (x 2 n + yl) - (x\ + yf) 
X = [x y] T . 


( 10 ) 


( 11 ) 


( 12 ) 


(13) 


(14) 


To estimate the receiver coordinates as X = [x y\ , the least 
square estimation method is employed to minimize the squared 


Euclidean distance as S = 


B — AX 


. Derivation of S 


A unique solution is finally calculated as 

X = (A t A)" 1 A t B. (16) 

IV. SIMULATION AND RESULT 

Three typical locations are selected to analyze the effect 
of multipath reflections. A (0, 0,1.2) represents a point at the 
corner of the room, where the scatterings and reflections are 
severe; B (4, 0, 1.2) represents a point at the edge of the room, 
right beside the wall, where reflections are medium; and C (4, 
4, 1.2) represents central point where the effect of multipath 
reflections becomes weak. Fig. 5 shows a birds-eye view of 
the system model. The circles show the locations of all the 
LED bulbs, and the three selected locations are marked with 
squares. We assume that a total number of 16 LED bulbs are 
installed on the ceiling of the room. The inner region is the 
area within the dashed lines in Fig. 5 while the rest of area is 
considered as the outer region. 

A. Impulse Response Analysis 

Considering the symmetrical property of the room, the im¬ 
pulse responses from the transmitter located at (3, 3, 3. 3) are 
investigated at the three selected locations. The contribution 
of the LOS and the first three reflections are shown in Figs. 
6 - 8 . 

Particularly, Fig. 6 demonstrates the impulse response of 
each reflection order at Location A. The impulse response 
amplitude of reflections is comparable to that of LOS incurring 
large positioning errors. Fig. 7 shows the impulse response 
of each reflection order at Location B. The amplitude of the 
reflections significantly decreases compared to Location A, 
and thus positioning accuracy is expected to be improved. 
As can be seen from Fig. 8, the LOS component almost 
dominates the total impulse response, and the amplitude of 
the reflections is negligible at Location C. Therefore, the 
positioning performance is expected to be less affected by 
multipath reflections. 
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Fig. 6. 


Fig. 7. 


Fig. 8. 
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x ^ q- 7 Impulse Response of Each Reflection Order at Location A 

—Contribution of LOS 
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—Contribution of 3-rd reflections 
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Impulse response of each reflection order at Location A. 



Time (unit: ns) 


Fig. 9. Received power from each LED bulb at Location A. 



Fig. 10. Received power from each LED bulb at Location B. 


Impulse response of each reflection order at Location B. 


xIO 


Impulse Response of Each Order of Reflections at Location C 


—Contribution of LOS 
— Contribution of 1-st reflections 
— Contribution of 2-nd reflections 
—Contribution of 3-rd reflections 
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Impulse response of each reflection order at Location C. 


B. Power Intensity Distribution Analysis 

As RSS information is applied to estimate the distance 
between a transmitter and receiver, the received power from 
each transmitter directly affects positioning performance. In 
this subsection, we investigate received power from different 
LED transmitters at each reflection order for the three selected 
locations. Figs. 9-11 present the highest six received power 
values for the selected locations inside the room in descending 
order. The received power of each reflection order at the 
corner point is shown in Fig. 9. As can be seen clearly, 
only for the first LED signal, the LOS power value is much 
greater than that of the reflections. However for the other 
five LED signals, the reflection components are comparable 
to the LOS component. The reflection components affect 
positioning accuracy since only direct power attenuation from 
the transmitter is considered in distance estimation. 


C. Analysis of Positioning Accuracy 

When distance between a transmitter and receiver is es¬ 
timated by the received power, the receiver coordinates are 
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Positioing Error(m) with Reflections 



Fig. 11. Received power from each LED bulb at Location C. 


Positioning Error without Reflections 



Fig. 12. Positioning error considering no reflections. 


Fig. 13. Positioning error considering reflections. 


1400 r 


Histogram of positioning error with no reflections 



0.005 0.01 0.015 0.02 0.025 

Positioning Error (m) 


0.03 


finally obtained by the positioning algorithm as stated in 
Section III. 

As a benchmark and in order to show the effect of multi- 
path reflections on the positioning accuracy, positioning error 
neglecting the reflected power is also calculated and shown in 
Fig. 12. As can be seen, the positioning error is low all over the 
room, and only slightly higher in the corner area. Fig. 13, on 
the other hand, shows the positioning performance considering 
the multipath reflections. It can be noted that each location 
of the room is affected by multipath reflections, especially 
the corner and edge areas. However, the positioning accuracy 
is satisfactory at the centre of the room as reflections are 
weak there. Fig. 14 and Fig. 15 present the histograms of 
positioning errors when neglecting and considering multipath 
reflections, respectively. When no reflections are considered, 
the errors only come from the thermal noise and shot noise 
which are small |[26l . In this case, most of the errors are within 
0.005 m which is an ideal scenario. However, reflections 
cannot be practically avoided, and thus they are a major 
concern in the positioning system impairing dramatically the 
system performance as shown in Fig. 15. In this case, most of 


Fig. 14. Histogram of positioning error considering no reflections. 


Histogram of Positioning Error with Reflections 

600|-T-T-T-T-T-T- 



Positioning Error(m) 

Fig. 15. Histogram of positioning error considering reflections. 
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TABLE III 

Positioning error neglecting/considering 

REFLECTIONS 


RMS error 

Reflections 
neglected (m) 

Reflections 
included (m) 

Loc. A 

0.0098 

1.6544 

Loc. B 

0.0019 

0.9966 

Loc. C 

0.0012 

0.1674 

Outer region (RMS) 

0.0059 

0.8173 

Inner region (RMS) 

0.0016 

0.2024 

Entire room (RMS) 

0.0040 

0.5589 


the positioning errors are below 1 m while at some locations, 
the error climbs up to 1.7 m. 

Table III compares the positioning error quantitatively when 
neglecting and considering the reflections. At Location A, the 
root mean square (RMS) error is 1.6544 m since the reflections 
are strong there. The effect of multipath reflections is medium 
at Location B while the positioning performance is the best 
at Location C. For the outer region, the RMS error is 0.8173 
m due to severe reflections while the RMS error of the inner 
region is 0.2024 m. The RMS error for the entire room is 
0.5589 m while it is only 0.0040 m when no reflections are 
considered. 

V. Discussion and Calibration Approaches 

As shown in Section IV, multipath reflections considerably 
affect the positioning accuracy, especially on the outer region. 
To improve system performance, three calibration approaches 
are proposed in this section. 

A. Nonlinear estimation 

Previous calculation is based on the linear least square 
estimation method as stated in Section III. Coordinates in © 
are obtained by looking for the optimized solution minimizing 
S. However, the mathematical deduction from © to © is not 
reversible. In other words, the optimized solution for © is not 
perfectly suitable for ®. Therefore, the linear least square 
estimation may introduce an error in positioning algorithm. 
Nonlinear least square estimation is proposed in this section 
as a calibration approach. Instead of finding the optimized least 
square solution for ©, solution for © is directly estimated 
by finding X = [x,y\ that minimizes 

S = y: {F(x, y, Xj, Di) - rf) , (17) 


where F(x , y ; Xi, yt) = y ( x- Xif + {y - yff . Employing 
the trust region reflective algorithm, an iterative process is 
applied to estimate X 11271 . Briefly, this algorithm works by 
first initializing an estimate as Xo, and then calculate the 
corresponding So- Second, several points surrounding Xo are 
substituted in ( fl7l >. and the one that minimizes S x is selected 
as Xt. After several iterative steps, receiver coordinates X 
will finally be obtained when S converges. In our approach, 
the value estimated with linear least square estimation in (fl4] i 


Positioing Error(m) with Reflections by Nonlinear Estimation 



Fig. 16. Positioning error with nonlinear estimation. 


Histogram of Positioning Error with Reflections by Nonlinear Estimation 



Fig. 17. Histogram of positioning error with nonlinear estimation. 

is selected as the initial value. Fig. 16 demonstrates the 
positioning error distribution with the nonlinear estimation 
approach. The histogram of the positioning errors is shown 
in Fig. 17. In this case, most of the errors are within 0.8 m, 
and only a few of them are over 1 m. The worst positioning 
error is just around 1.5 m. 

Table IV compares the performance of the linear and non¬ 
linear estimation methods. As can be seen, the nonlinear 
estimation outperforms its linear estimation counterpart. Par¬ 
ticularly, for the outer region where the reflections are severe, 
the nonlinear estimation reduces the RMS error to 0.6871 m. 
The RMS error also decreases to 0.1401 m and 0.4642 for the 
inner region and the entire room, respectively. 

B. Selection of LED Signals 

The received power decreases when distance between a 
transmitter and receiver increases. As shown in Fig. 5-7, the 
reflections contribute more to the total received power when 
the signal is from a farther LED transmitter and brings a larger 
error in the distance estimation. 






















TABLE IV 

RMS ERROR WITH LINEAR/NONLINEAR ESTIMATION. 


RMS error 

Outer 

Region(m) 

Inner 

Region(m) 

Entire 
Room (m) 

Linear 

estimation 

0.8173 

0.2024 

0.5589 

Nonlinear 

estimation 

0.6871 

0.1401 

0.4642 


TABLE V 

RMS ERROR WITH TRANSMITTER SELECTION APPROACH. 


RMS error 

Outer 

Region(m) 

Inner 

Region(m) 

Entire 
Room (m) 

6 LEDs (Linear 
estimation) 

0.6760 

0.1640 

0.4616 

6 LEDs (Nonlinear 
estimation) 

0.6016 

0.1112 

0.4046 

5 LEDs (Linear 
estimation) 

0.5472 

0.1251 

0.3722 

5 LEDs (Nonlinear 
estimation) 

0.5272 

0.0851 

0.3527 

4 LEDs (Linear 
estimation) 

0.4838 

0.0933 

0.3259 

4 LEDs (Nonlinear 
estimation) 

0.4828 

0.0849 

0.3240 


In this subsection, a signal selection approach is proposed, 
i.e., the receiver only selects strong signals for the coordinate 
estimation. In the numerical analysis, the six, five and four 
strongest LED signals are selected, and the results are shown 
in Table V. By removing the signals affected considerably by 
multipath reflections, the positioning accuracy is improved. 
The total RMS error decreases to 0.4046 m, 0.3527 m and 
0.3240 m respectively for the cases when the six, five and 
four strongest LED signals are selected. Specially, when the 
four strongest LED signals are selected, the RMS error is 
0.4828 m for the outer region and 0.0849 m for the inner 
region. As can be noted, the outer region is improved more 
than the inner region. For the sake of presentation, only the 
best results are shown in Figs. 18 and 19. Fig. 18 shows the 
positioning error distribution when the four strongest LED 
signals are selected for distance calculation and the nonlinear 
estimation is applied to obtain receiver coordinates. Fig. 19 
presents the corresponding histogram of positioning errors. It 
can be seen from Fig. 19 that many of the locations have 
positioning error which is less than 0.4 m, and only a few 
locations have positioning error that is larger than 0.8 m. 

C. Deceasing the Distance between LED Bulbs 

When LED bulbs are installed in a dense layout (i.e., 
the distance between LED bulbs is reduced, and a greater 
number of LED bulbs are used), the light intensity distribution 
becomes more uniform for the entire room, and therefore, the 
positioning accuracy is improved. 

Table VI shows the RMS error where distance between 
the LED bulbs decreases to 1.5 m, and 25 LED bulbs are 
installed in total. With no LED signal selection, the entire 
RMS error is 0.3121 for the nonlinear estimation. The RMS 
error decreases to 0.2922 m and 0.2699 m when six and five 
LED signals are selected, respectively. However, when the four 
strongest LED signals are selected, the receiver coordinates 


Positioing Error(m) with 4 LEDs by Nonlinear Estimation 

1.2 



Fig. 18. Positioning error (4 selected LED signals and nonlinear estimation). 


Histogram of Positioning Error with 4 LEDs by Nonlinear Estimation 



Fig. 19. Histogram of positioning error (4 selected LED signals and nonlinear 
estimation). 


cannot be obtained by the linear estimation in some locations. 
The shorter distance between the LED bulbs increases the 
probability that three of the selected LED bulbs are in the 
same row. In this scenario, there is a singularity in Matrix A, 
and the coordinates estimation fails. If this case happens, the 
average of the detected LED bulb coordinates will be applied 
as the initial value instead of the linear estimated value. As 
this initial value is not as accurate as the linear estimated 
one, the positioning error increases as shown in Table VI. 

Figs. 20 and 21 show the positioning performance for the 
best scenario, i.e., 5 LED signals selected for the nonlinear 
estimation. As can be seen from Fig. 20, for the most of the 
inner area, the positioning performance is satisfactory while 
at the edges and corners of the room, there are some locations 
with large positioning errors. From Fig. 21, it can be observed 
that most of the positioning errors are within the range of 0.5 
m demonstrating improvement in the positioning performance 
and usefulness of the proposed method. 
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TABLE VI 

RMS ERROR WITH 1.5 M DISTANCE BETWEEN THE LED 
BULBS. 


RMS error 

Outer 

Region(m) 

Inner 

Region(m) 

Entire 
Room (m) 

All LEDs (Linear 
estimation) 

0.6049 

0.1093 

0.4064 

All LEDs (Nonlinear 
estimation) 

0.5310 

0.0861 

0.3553 

6 LEDs (Linear 
estimation) 

0.4400 

0.0916 

0.2976 

6 LEDs (Nonlinear 
estimation) 

0.4362 

0.0729 

0.2922 

5 LEDs (Linear 
estimation) 

0.4106 

0.0691 

0.2751 

5 LEDs (Nonlinear 
estimation) 

0.4014 

0.0742 

0.2699 

4 LEDs (Linear 
estimation) 

NA 

0.0571 

NA 

4 LEDs (Nonlinear 
estimation) 

0.4251 

0.0729 

0.2850 


Positioing Error(m) with 5 LEDs by Nonlinear Estimation 



Y ( m > 0 “ X(m) 


Fig. 20. Positioning error with 1.5 m distance between LED Blubs (5 selected 
LED Signals and Nonlinear Estimation). 


Histogram of Positioning Error with 5 LEDs by Nonlinear Estimation 



Fig. 21. Histogram of Positioning Error with 1.5 m Distance between LED 
Blubs (5 selected LED Signals and Nonlinear Estimation). 


In this paper, an indoor visible light positioning system 
taking account of multipath reflections has been investigated 
for a typical room where the impulse response is obtained 
employing CDMMC approach. In particular, the impulse re¬ 
sponse of each order has been calculated for three locations 
representing the corner, edge and center points. The received 
signal power distribution at three locations with different de¬ 
grees of multipath reflections has been analyzed. Furthermore, 
the positioning error in the entire room has been calculated 
and compared with the previous works where reflections are 
not taken into account. Comparison has shown that multipath 
reflections considerably decrease the positioning accuracy, 
especially for the outer region. Three calibration approaches 
have been proposed to improve the positioning performance. 
Particularly, we have shown that by employing nonlinear 
estimation, the positioning error is decreased. A selection on 
the LED signals has been also proposed to remove the signals 
that are severely affected by multipath reflections. Finally, 
we have shown that a dense layout of LED bulbs results in 
a more uniform light intensity distribution, and therefore a 
higher positioning accuracy can be achieved. 
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